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Abstract. Supermassive black holes inhabit galactic nuclei, and their presence influences 
in crucial ways the evolution of the stellar distribution. The low-density cores observed in 
bright galaxies are probably a result of black hole infall, while steep density cusps like those 
at the Galactic center are a result of energy exchange between stars moving in the gravita- 
tional field of the single black hole. Loss-cone dynamics are substantially more complex in 
galactic nuclei than in coUisionally-relaxed systems like globular clusters due to the wider 
variety of possible geometries and orbital populations. The rate of star-black hole interac- 
tions has begun to be constrained through observations of energetic events associated with 
stellar tidal disruptions. 



1. Introduction 

The association of supermassive black holes 
(SBHs) with galactic nuclei began even be- 
fore their observation a l confi rmation, when 
IZel'dovich & Novikovl ill 964 and Sal peteJ 
ill 964 1) first discussed the growth of mas- 
sive objects at the centers of galaxies. Since 
then, the mutual interaction of SBHs and 
stars has played a central role both in the 
detection and mass d etermination of SBHs 
llFerrarese & Ford 2005). and also in our the- 
oretical understanding of how nuclei form and 
evolve. The realization that galactic spheroids 
grow through mergers added an interesting 
complication to this picture, since galaxy 
mergers imply the formation o f binary SBHs 
llBe gelman. Blandford & ReesI f!98Q.) . which 
can influence the stellar distribution on larger 
spatial scales than single SBHs, and which may 
be direc tly or indirectly c onnected with nuclear 
activity llKomossal2003h . 



2. Binary Black Holes and Cores 

Galaxy mergers bring SBHs together, and 
unless the binary SBHs themselves coa- 
lesce, ejections will occur when a third 
SBH is deposited into a nucleus containing 
an uncoalesced binary ( Mikkola & Valtoneii| 
.1990 : Volonteri, Haardt & Madau 2003). Such 
events can not be too frequent or the tight, em- 
pirical correlations between SBH masses and 
galaxy properties ( Ferrarese & Merritt 2000; 
iGraham etaf] l2(K)ll: Marconi & Hunt 2003) 
would be violated ( Haehnelt & Kauffmanril 
I200I . 

Ultimately, coalescence is driven by emis- 
sion of gravitational waves, but this process 
only becomes efficient when separations fall 
below ~ 10"^ pc. By comparison, the natural 
separation between two SBHs of mass nii and 
OT2 in a galactic nucleus is 
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Fig. 1. Short-term evolution of binary black holes 
in spherical galaxies. Upper panel: Thick (black) 
line shows infall of a hole of mass mi into a nu- 
cleus containing a hole of mass m\ = 10m2. Thin 
(red) line is the orbital decay predicted by the dy- 
namical friction equation, assuming an unchanging 
galaxy. The vertical dashed line indicates the time 
when a x a/,. Lower panel: Evolution of the mass 
deficit, i.e. the mass displaced by the binary, in units 
of the combined mass of the two black holes. The 
evolution slows at a w a;, due to ejection of stars 
on orbits that intersect the binary, and thereafter the 
evolution rate is strongly dependent on the number 
A' of particles used to represent the galaxy (here, 
N = 2x 10^). For much larger A', the binary would 
stall at a » fl;, (see Fig. 2). 



where M, - m\ + m2, i-i - m\m2lM, is the 
reduced mass, q - m2lmi < 1 is the mass ra- 
tio, M..8 - M.I Iff Mq, and 0-200 is the nuclear 
velocity dispersion in units of 200 km s"'. A 
binary with semi-major axis a ^ ah'is "hard," 
i.e. its binding energy per unit mass is ~ cr^ 
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Fig. 2. Long-term evolution of binary black holes 
in spherical galaxies. Each curve shows the evolu- 
tion of the inverse semi-major axis of an equal-mass 
binary in a spherical galaxy with an initially power- 
law nuclear density profile. Curves were obtained 
by averaging multiple A'-body integrations with the 
same A'. As A' increases, the binary's hardening rate 
drops, implying very slow evolution past a « a;, in a 
real galaxy (adapted from Merritt, Mikkola & Szell 
2006). 



and its total binding energy £/, is a fraction 
~ M./Mgai ~ lO"-' that of the host galaxy. 

The separation a/, is natural one since a bi- 
nary with a < ai, efficiently ejects stars on 
intersecting or bits (Mikkola & Valtonea 1992i: 
IOuinlanlll996h . In a spherical or axisymmetric 
galaxy, the total mass associated with stars on 
such orbits is small, < M,, and once the binary 
has ejected these stars, its evolution greatly 
slows. This is the "final parsec problem." In A^- 
body simulations (e.g. Milosavljevic & Merritt 
2001, Makino & Funato 2004), the binary of- 
ten does evolve past a ^ o/, but at a rate that 
is strongly A^-dependent (Figs. I1I2> . This re- 
flects the fact that processes like gravitational 
scattering and Brownian motion that can refill 
a binary's loss cone are m^, -dependent. 

Since the efficiency of loss-cone repopula- 
tion in real galaxies is uncertain and may be 
low, it is interesting to ask what the observable 
consequences would be of binaries stalling at 
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Fig. 3. Observational determination of mass 
deficits. Plot shows the observed surface brightness 
profile of NGC 3348. The dashed line is a Sersic 
model fit to the large-radius data. Solid line is the 
fit of an alternative model, the "core-Sersic" model, 
which fits both the inner and outer data well. The 
mass deficit is illustrated by the area designated 
"depleted zone" and the corresponding mass is 
roughly 3 x IO'^Mq (from Graham 2004). 

a ^ ah. Using the M, - cr relation, Eq. Q can 
be written 

a/i,pc ~ 0.3 ^0.1 o"200 (2) 

with q'o.i = q/Q.l. "Dual" SBHs are observed 
in a hand ful of interactin g gala xies and binary 
quasars (Mortlock et al. 1999; Komossa et al.l 
[2003:.Ballo et al..2004) . but always with much 
larger separations, > 1 kpc. True binaries with 
a * ai, would be difficult to resolve outside 
the Local Group even if both SBHs were lu- 
minous. A number of active galaxies show evi- 
dence for variability with periods similar to bi- 
nary orbital periods, but the variability may be 
a scaled-up version of the quasi-periodic oscil- 
lations associated with stellar-mass blac k holes 
and n ot the signature of a binary SBH llUttlevI 
120051) . 

An indirect way to track the evolution of a 
binary SBH is via its effects on the structure 
of the nucleus. The infalling hole begins dis- 
placing stars when it first becomes bound to 
the second, at a separation ~ r/, = GM,/cr^ ^ 
lOM.^go-jQo pc X ISM'-^ f pc. Evolution from 



rh to Oh is rapid unless the mass ratio is ex- 
treme, and the "damage" done by the binary 
to the nucleus can be robustly estimated from 
A^-body simulations, independent of uncertain- 
ties about the efficiency of loss-cone refilling 
(Fig. 0. A standard definition of the damage 
is the "mass deficit" M^ief, the difference in 
integrated mass between the observed density 
pro file and the p rimordial (pre-merger) pro- 
file (iMilosavhevic et al 2002) . A^-body simu- 
lations like those in Fig. 1 show that Mdefia - 
a/,) ~ 0.5(mi + 1712), decreasing only weakly 
with 1112, and nearly independent of the form of 
the initial density profile. This is a consequence 
of the scaling of Ei, with mi -i- 1112 discussed 
above. An interesting corollary is that the mass 
deficit after N mergers should be ~ 0.5M,N, 
with M, the final (i.e. current) mass of the 
SBH, almost independent of the details of the 
merger tree. 

Bright (Mb < -20.5) elliptical galaxies 
always exhibit cores, regions near the cen- 
ter where the luminosity profile is nearly flat 
(Ferrarese et al. 2006 and references therein). 
Estimating Mjp/ in a real galaxy is not straight- 
forward however, since one does not know 
what the pre-existing stellar distribution was. A 
conservative approach is to fit a smooth func- 
tion, e.g. Sersic's law, to the outer luminosity 
profile and extend it inward; the mass deficit 
is then defined in terms of the differential pro- 
file (Fig.O. Mass deficits derived in this way 
are found to be roughly equal to M, ( GrahanJ 
12004: .Ferrare se et al. 2006) This is consistent 
with the relation given above if AT ~ 2, 
in reasonable accord with hierarchical galaxy 
formation models. Alternatively, if the binary 
evolves all the way to coalescence and if gravi- 
tational waves are emitted anisotropically dur- 
ing the final plunge (the "rocket" effect; Favata 
et al. 2004), the displaced hole will transfer en- 
ergy to the nucleus before falling back to the 
center, increasing the mass deficit b y some ap- 
preciable fraction of M, ( Merritt et'ani2004l 
Boylan-Kolchin et al. 2004- 

Evolution past a ^ an may or may not 
leave an additional mark on a nucleus. In small 
dense galaxies like M32, two-body (star-star) 
relaxation can repopulate a binary's loss 
cone at interesting rates, possibly driving the 
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Fig. 4. Efficient merger of binary SBHs in barred galaxies. The left panel shows the evolution of the bar 
amplitude in a set of initially axisymmetric galaxy models with different degreess of rotational support 
(measured by the parameter too)- For a>o > 0.6 the models are bar-unstable. The right panel shows the 
hardening rate of a massive binary at the center of the galaxy. In the barred (triaxial) models, the binary's 
hardening rate is independent of particle number A', while in the axisymmetric models it falls roughly as 
yV"' indicating a coUisionally-resupplied loss cone (adapted from Berczik et al. 2006) 



binary to coalescence in < 10 Gyr ( lYdl2002t 
iMilosavlievic & Merri3 '2Q03'). A^-body simu- 
lations like those in Fig. |2] show that the mass 
deficit increases as M^/ef oc ln(fl/,/fl) in such 
cases, so that the observed value of M^ef does 
not contain much information about the final 
value of a. The bright elliptical galaxies that 
are observed to contain cores have relaxation 
times that are much too long for coUisional loss 
cone repopulation to occur however A more 
efficient way to drive a binary to coalescence 
is to imbed it in a strongly n on-axisymmetric 
(barre d or triaxial) galaxv ([Memtt^PooE 
|2224t iHoUev-Bockel mann & Sigurdssori 
|200d). Triaxial potentials can support large 
populations of "centrophilic" orbits, allowing 
a binary to harden well past a/, by interacting 
with stars even in the absence of collisional 
loss-cone refilling (Fig. |3. Most of the stars 
supplying mass to the binary under these 
circumstances have orbital apocenters » r/, 
and their depletion has almost no effect on the 
nuclear density profile. 

A binary with a <K a/, ejects stars with 
velocities large enough to expel them from 
the galaxy. A handful of hypervelocity stars 
(HVSs) have been found in the Galactic halo 



that are candidates for e j ection (lBro\yn et alJ 
120051 iHirsch etal J 1200.4 iBrown et aU MOOfh. 
These stars must have left the nucleus < 10** 
yr ago however, implying that a binary SBH 
was present until quite recently. Other models 
can produce HVSs in the absence of a second 
SBH (e.g. Hills 1988; Gualandris et al. 2005). 
Ejection by binary SBHs might also be re- 
sponsible for populations like the intergalactic 
planetary nebulae observed in the Virgo clus- 
ter ( Hollev-Bockelmann et al. 2005). However, 
such models require that the binary SBHs 
ejecting the nebulae harden well beyond a = 
fl/,, and also that the hardening is driven by 
star-binary interactions; both assumptions are 
open to question, particularly in the case of the 
brightest Virgo galaxies. 

The nuclei of merging galaxies are ex- 
pected to contain some of the largest concen- 
trations of dense gas in the universe. Binary- 
gas simulations fall into two categories de- 
pending on whether the gas is modelled as 
"hot," i.e. having a specific kinetic energy com- 
parable with that of the stars, or "cold," e.g. 
in a thin disk. The first case is much easier to 
treat; if the mass in gas within ~ a/, is compa- 
rable to M„ the additional component of dy- 
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namical friction drag can bring t he two SBHs 
together in m uch less than 10^ yr llEscala et al.l 
l200il2005l) . However it is not clear that such 
massive, dense accumulations of hot gas can 
be sustained. In the case of binary-disk inter- 
actions the difficulties are computational, e.g. 
resolving the gas on scales smaller than the ra- 
dial wavelength of the density waves induced 
by the binary. Much progress in this area is to 
be expected in coming years. 

3. Nuclear Equilibria 

Once the massive binary has coalesced (or oth- 
erwise stopped interacting with stars), the stel- 
lar distribution in the nucleus can evolve to a 
steady state. An important distinction can be 
made between colUsional nuclei, which have 
two-body (star-star) relaxation times less than 
10 Gyr, and collisionless nuclei, with Tr > 10 
Gyr. Fig. |5] shows that Trifh) drops below 10 
Gyr in spheroids roughly as faint as the bulge 
of the Milky Way. 

The morphology of collisionless nuclei is 
constrained only by Jean's theorem, i.e. by the 
requirement that the density be representable 
as a superposition of orbits integrated in the 
self-consistent potential. Collisionless equilib- 
ria can take many form s, including triaxial nu - 
clei with chaotic orbits llPoon & Merriti2004 . 
as well as the simpler, axisymmetric models 
that are the basis for mo st SBH mass estima- 
tion llGebhardt et alJl2003h . A recent review is 
given by Merritt (2006). 

In a collisional nucleus, on the other hand, 
exchange of energy between stars drives the 
phase-space density toward a characteristic 
form. Here we discuss the more limited class of 
solutions associated with collisionally-relaxed 
nuclei. 

Gravitational encounters drive the velocity 
distribution of stars around a black hole to- 
ward a Maxwellian on a time scale of ~ T,., 
but a Maxwellian velocity distribution implies 
an exponentially divergent mass near the hole. 
The existence of a region close to the black 
hole in which stars are captured or destroyed 
prevents the nucleus from reaching precise 
thermal equilibrium. The density must drop to 
zero on orbits that intersect the hole's event 
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Fig. 5. Nuclear relaxation times in a sample of 
early-type galaxies and bulges. The size of the sym- 
bols reflects the degree to which the SBH's influence 
radius was resolved; filled symbols have O^i^ > 9abs 
(resolved) and open circles have < 0o4, (unre- 
solved) (from Merritt & Szell 2005).' 

horizon or that pass within the tidal disruption 
sphere at r = r,; the latter radius is most rele- 
vant here since galaxies with collisional nuclei 
probably always have M, < IO^Mq (Fig. |5} 
and hence r, > IGM.jc^. 

For a single-mass population of stars mov- 
ing in the Keplerian potential of a black hole, 
the steady-state solution to the isotropic, orbit- 
averaged Fokker-Planck equation in the pres- 
ence of a central "sink" is approximately 

f{E) = ME\'^\ p(r)^por-''\ (3) 
\Eh\ « |£| « \E,\, r, « r « r,, (4) 

jRahcall & Wolfl Il976t iTjghtman & ShanirnI 
Il977h . Here~ = -GM./r, and Eh = 
-GM.jrh- This is a "zero-flux" solution, i.e. it 
implies = at \Eh\ « \E\ « \Et\ where Fe 
is the encounter-driven flux of stars in energy 
space. The actual flux is smaU but non-zero, of 
order 

Ff a oc n. (5) 

Trin) 

In other words, the flux is limited by the rate 
at which stars can diffuse into the disruption 
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sphere at r,. This flux is "small" in the sense 
that the one-way flux of stars in or out through 
a surface at r, «; r <s; r;, is much greater; ex- 
cept near r,, the inward and outward fluxes al- 
most cancel. 

The Bahcall-Wolf solution has been 
verified in a nu mber of studies ba sed o n 
Fokker- Planck llCohn & KulsrudI Il978h . 
fluid ( Amaro-Seoane et al. '20041), or 
Monte-C ai-lo ( Mai'chant & Shapiro 11980^ 
[ Duncan & Shapiro 1983; Freitag & Benz 
l2002h approximations. Most recently, ad- 
vance s in computer hardware ( Makino e t alJ 
l2003h have made it possible to test the 
Bahcall-Wolf solution via direct A^-body inte- 
grations, avoiding the appr oximations of the 
Fokker-Planck formalism (Pretoetal. 2004; 
[Baumgardt et al. 2004; Merritt & Szell 2005). 
Another advantage of this approach is that it 
can easily deal with complex initial conditions 
and nonspherical geometries, such as those set 
up by a binary SBH. 

Fig-Elshows a simulation of cusp destruc- 
tion (by a binary SBH) followed by cusp re- 
generation (via the Bahcall-Wolf mechanism). 
The Bahcall-Wolf cusp rises above the core in- 
side a radius ~ 0.2r/, and reaches its steady- 
state form in a time 2 - 3r,(0.2r/,). In the Milky 
Way nucleus, 0.2r/, » 0.7 pc and 7, (0. 2r/,) » 6 
Gyr (assuming = 0.7 Mq). In fact, as Fig.|6l 
shows, the stellar distribution at the Galactic 
center is close to the Bahcall-Wolf form. 

Figs. |5] and |6l suggest that Bahcall-Wolf 
cusps should be present in many galaxies with 
luminosities similar to that of the Milky Way 
bulge, or fainter, assuming that they contain 
SBHs. M32 fits these criteria; th e mass of its 
SBH is uncertain by a factor ~ 5 dValluri et al.l 
I2004h but even if the largest possible value 
is assumed, 0.2r/, w 0.8 pc w 0.3" would 
only barely be resolved. NGC 205 clearly con- 



4. Loss-Cone Dynamics 

As discussed above, the existence of a region 
r < r, close to the SBH where stars are de- 
stroyed has a profound influence on the steady- 
state distribution of stars even at r » r,. 
Loss of stars is also important because of its 
observational consequences: tidally disrupted 
stars are expected to produce X- and UV ra- 

diation with luminosities of ~ K)^"* erg s~' 

( iMeszaros & Silk"1977; Rees 1990). 

Once a central SBH has removed stars 
on orbits that intersect the disruption sphere 
r = r,, continued supply of stars to the 
SBH requires some mechanism for loss- 
cone repopulation. The most widely dis- 
cussed mechanism is gravitational encounters, 
which drive a diffusion in energy and angu- 
lar momentum. The latter dominates the loss 
rate, since even a small percentage change 
in angular momentum can put a typical 
star into the loss cone (' Frank&ReesllI976l 
lUghtman & Shapiro 1977). 

A large fraction of the stars within the 
SBH's influence radius will be deflected into 
r, in one relaxation time, i.e. the loss rate is 
roughly {M,/mi,)/Tr{ri,). In a collisional nu- 
cleus with M. « lO^'Mo, this is ~ 10^(10'" 
yr)^ lO"'* yr-i. A slightly more careful calcu- 
lation for the tidal disruption rate in a p oc r^^ 
nucleus gives 



N ^ Ix lO^V"' 



\70kms" 



7/2 



/ M. ^ 



-1/3 



(6) 



where and are t he mass and radius of th e 
tidally disrupted stars dWang & Merrittll2004h . 
Assuming p ~ r^^ is reasonable since this 
is approximately the slope observed with the 
Hubble Space Telescope at r > r/, in the low- 
lumi nosity galaxies that dominate the event 
rate rTehhardtetall1996h. 

A handful of X-ray flares have been de- 



no dvnamical evidence of a massive black dKomossa & Greinei 




Il999t 


hole ( Valluri et al. 


2005). AU other ealaxies in l2004l iHaloern et al. 




2004ir 



which r;, is resolved have nuclear relaxation 
times » 10Gyr(Fig.|5j. 



iKomossa et al.1 
The mean dis- 
ruption rate computed from this small set of 
events is very uncertain but is probably con- 
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Fig. 6. Cusp regeneration. Curves show the evolution of the spatial (a) and projected (b) density profiles 
of an A'-body model in which an initial, p ~ r"'-' density cusp (dotted line) was destroyed by infall of a 
second black hole of mass m2 = 0.5mi. Blue (lower) line shows the density after infall, and green, red and 
black lines show the evolving density after the two black holes have been merged into one; the final time 
is ~ 10 Gyr if scaled to the Milky Way. Symbols in (b) show the observed surface density of stars near the 
Galactic Center (Genzel et al. 2003). Dashed lines have logarithmic slopes of -1.75 (a) and -0.75 (b), the 
Bahcall-Wolf "zero flux" solution (from Merritt & Szell 2005). 



sistent with Eq. ^ llDonlev et alJl20Q2h . The 

nuclei in which tidal disruptions have occurred 
appear to remain luminous for 1 - 10 yr after 
disruption and possibly longer; in three of 
the observed events, the luminosity decay ap- 
proximately obeys the oc r^l^ dependence 
predicted if emission is produced during the 
fallback of stellar debris onto an accretion disk 
llPhinnevll989HEvans & Kochanekll98l . 

C lassical loss cone theory 
JUghtman & Shapiro 1 9771; ICohn & Kulsrud 
Il978h . from which expressions like Eq. (|6j 
were derived, was directed toward understand- 
ing the observable consequences of massive 
black holes at the centers of globular clusters, 
which are many relaxation times old. As Fig.|5] 
shows, few galactic nuclei are expected to be 
much older than T,-, and loss cone dynamics in 
galactic nuclei can therefore be very different 
than in globular clusters. For instance, in a 
nucleus that until recently contained a binary 
SBH, orbits of stars with pericenters rperi < fl/i 
will have beeen depleted. The time required 
for gravitational encounters to repopulate 



these orbits is ~ {aiJri,)T,- ^ {m2lmi)T,-. For 
m2/mi = 0.1, this time exceeds 10"^ yr in most 
galax ies with M. > 10** Mq ( Merritt & WansJ 
120051) . Until the phase-space gap is refilled, 
tidal disruption rates can be much lower than 
in a collisionally-relaxed nucleus. 

On the other hand, as Fig.l^implies, loss- 
cone repopulation in a non-axisymmetric (tri- 
axial or barred) nucleus can be much more ef- 
ficient than repopulation due to gravitational 
encounters in the spherical geometry, due to 
the presence of centrophilic (box or chaotic) 
orbits in the triaxial geometry. Numerical in- 
tegrations of centrophilic orbits show that the 
rate at which a single star experiences near- 
center passages wit h pericenter distanc es < d 
is proportional to d llGerhard & Binnevl fr985). 
In a p ~ r^^ nucleus, the implied rate of sup- 
ply of stars to the event horizon of the SBH is 
~ fccr' IGc^ where fc is the fraction of orbits 
that are centrophilic. Even for modest values 
of fc {~ 0.1), this collisionless mechanism can 
supply stars to the SBH at higher rates than col- 
lisional loss-cone repopulation, particularly in 
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galaxies with M, > 10^ Mq in which two-body 
relaxa tion times are very long llMerritt & PoonI 
l2004ft . In fact, loss rates in the triaxial geome- 
try can approach the so-called "full loss cone" 
feeding rates in spherical galaxies, which were 
invoked, in an early model, to explain QSO 
activity jHills" '1975', 'Young et all Il977l) and 
which hav e recently been revived ( Zhao et al. 
120021 iMiralda-Escude & KoUme ier 2005). 

These arguments suggest that the mean rate 
of stellar tidal disruptions in galactic nuclei is 
poorly constrained by pure theory. One way 
to observationally constrain the disruption rate 
is via the X- ray luminosity fu n ction of ac- 
tive galaxies (Ueda et al.' l 2003t iBarger et alJ 
[2005; Hasinger et al. 2005]). Assuming the 
Lx r^^^ time dependence discussed above 
for individual events, and convolving Eq. (|6} 
with the SBH mass function, one concludes 
llMilosavhevic. Merritt & Hoi 1200 6') that tidal 
disruptions can account for the majority of X- 
ray selected AGN with soft X-ray luminosities 
below - 10"*^ - 10'*'*ergs-i. 

Nearer to home, it might be possible to 
search for "afterglows" of the most recent tidal 
disruption event at the Galactic center, which 
could plausibly have occurred as little as ~ 
10^ yr ago (Eq.|6j. Possible examples of such 
signatures include X-ray flourescence of giant 
molecular clouds (Sunyaev & Churazov 1998) 
and chang es in the surface prop erties of irradi- 
ated stars iTimenez et alJ200^ . 
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